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We derive the special and general relativistic hydro-dynamic equations of motion for 
ideal fluids from a variational principle. Our approach allows to find approximate 
solutions, whenever physically motivated trial functions can be used. Illustrating 
this, a Rayleigh-Plesset type equation for the relativistic motion of a spherical 
(QGP) droplet is obtained. Also the corresponding general relativistic effective 
Lagrangian for spherically symmetric systems is presented. 
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1 Introduction 

First applications of relativistic hydrodynamics to the process of multiparticle 
production in high-energy, hadronic collisions were studied by Fermi and Lan- 
dau in the early 1950'sBtl Recently, the relativistic motion of fluids has been 
studied, extensively with respect to the analysis of relativistic heavy-ion col- 
lisions, a The hydrodynamic description of high-energy hadronic and nuclear 
collisions has been successful in reproducing many characteristic features of 
these processes, such as multiplicity and transverse energy distributions. 

Theoretically, however, the foundation of the hydrodynamical picture for 
these processes is not fully understood. The application of hydrodynamics im- 
plicitly assumes local thermal equilibrium via an equation of state. Therefore, 
the relaxation time and the mean free path should be much smaller compared 
to, respectively, the hydrodynamical time scale and spatial size of the system, 
which may be satisfied best in collisions of heavy nuclei. 
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From a kinematical point of view, besides the equation of state, the equa- 
tions of hydrodynamics simply express the conservation laws of energy, mo- 
mentum, and possibly various charges. Thus, for processes where flow is an 
important factor, a hydrodynamic description seems natural to begin with. 
Effects of finite relaxation time and mean-free path may be implemented later 
on, viscosity and heat conductivity in particular, or^ome simplified transport 
equations might be applicable, for example, see Ref.Q. 

Rclativistic hydrodynamics is a local description of the conservation laws, 
i.e. in terms of the energy-momentum tensor, 

= , (1) 

where we momentarily assume cartesian coordinates. It is difficult to solve this 
set of coupled partial differential equations in any generality. Few analytical 
solutions are known, but even for simple geometries, like one-dimensional or 
spherically symmetric cases, one often has to resort to numerical solutions. 
Under these circumstances a physically motivated and more qualitative ap- 
proach, such as via a variational formulation, should provide useful insight, 
before or without full scale computer simulations. 

Another important arena for relativistic hydrodynamics is found in cosmol- 
ogy and high-energy astrophysics, such as the gravitational collapse of stellar 
cores forming a neutron stator a black hole, relativistic blast waves in models 
of gamma ray bursts, etcfl'tl The assumption of local thermodyrLamical equi- 
librium is considered to be well justified in most cases, cf. Ref.Ll. However, 
often the fluid dynamics and the local changes of the gravitational field have to 
be considered simultaneously. Again, this is generally very complicated and, 
therefore, a variational scheme seems extremely helpful. 

Frequently even the equation of state of matter is not known precisely. 
Rather, hydrodynamics is employed to infer its properties involved in the pro- 
cess. In such cases we need to describe the flow, which characterizes the dy- 
namics, assuming a specific equation of state. 

In conclusion, a variational method is highly desirable which allows to 
solve the dynamical equations of the system more effectively. Such an effective 
Lagrangian and variational principle has previously been employecLto study 
the effect of local turbulent motion on the supernova mechanism.ErEI It has 
also been useful in deriving the Rayleigh-Plesset equation and various gener- 
alizations for the gas bubble dynamics in a liquid.clilj Furthermore, the global 
features of high-energy hadxonic and nuclear collisions are described well by 
the fireball model, see RefB for a review. We aim at a dynamical scheme 
which improves the simplest fireball model in the direction of a realistic hydro- 
dynamical description, incorporating the QCD phase transition and important 
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dissipative processes. 

In the present work, we reformulate Special and General i?clativistic (SR 
and GR, respectively) hydrodynamics as a variational principle. We derive 
a SR generalization of the Raylcigh-Plesset equation in this approach. Fur- 
thermore, we obtain the GR effective Lagrangian for spherically symmetric 
systems. 



2 The Variational Approach 

2.1 Variables and Constraints 

While the variational formulation of nonrelativistic hydrodynamics has been 
established for quite a while, we consider SR and GR equations of motion 
here. The most natural variable describing the flow is the velocity field, v(f, t). 
Introducing the four- vector u^(x), we define instead the variational variables, 

m° = 7 , u = 7# , (2) 

i.e. functions of r and t, subject to the constraint, 

«/ = ! • (3) 

The flow of matter provokes local changes in the occupied volume. Considering 
the case of a conserved chargelike quantity, say the baryon number, let its local 
density in the comoving frame be n. Then, 

(nu^) = • (4) 

The energy of the matter is given by: 

E = sV , V = — , (5) 
n 

where e denotes the energy density and V the local specific volume. Assuming 
local equilibrium, we have the thermodynamical relations, 



{w) 8 = - ^ (l) s = £ ?' (6) 



where S is the entropy and p the local pressure. 
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2.2 The SR Action and Hydrodynamic Equations of Motion 



Let us first consider a single particle of rest mass mo with the action, 

Ip = —mo J ds = — J d^7 _1 mo = — J dt edV , (7) 

where 7 _1 = (1 — i; 2 ) 1 / 2 ; here we assumed an infinitesimal rest frame volume 
dVo of the particle, such that mo = edVo — ejdV in the laboratory frame, 
where the particle has velocity v. 

This consideration leads to the action for a fluid, Im = — j d^x e, which 
is considered here as an aggregation of infinitesimal volume elements. 

Varying the action with respect to the four- velocity field, ii M — > u M + Su^, 
and the density distribution, n — > n + Sn, we have to respect the constraints, 
Eqs. (Q), Introducing Lagrange multipliers, we obtain the fluid action, 

I M = Jd?x{-e(n)+£(x)dr[nu»] + l{(x)[u»ur-l]} , (8) 

and the variational principle is stated as: SIm — 0, for arbitrary variations of 
n, f, C- 

In practical applications of the variational principle it will be convenient 
to parametrize u M , n such that the constraints are automatically satisfied and 
the Lagrange multipliers £, £ eliminated from the beginning. 

We remark that the following derivation of the hydrodynamic equations 
of motion is equally valid in general coordinate systems. Then, the partial 
derivative 9 M should be replaced by the appropriate covariant derivative in the 
action and correspondingly in the following. Furthermore, the volume element 
d 4 x has to be replaced by the invariant volume element y / —gd 4 x$3 

Performing first a partial integration of the second term of the action, 
Eq. m), and then the four variations, we obtain the set of equations: 



u 2 = 1 , (9) 

a^nu") = , (10) 

|£+^ e = 0i (11) 

nd^ -(u^O , (12) 



from which the Lagrange multipliers can be eliminated. 

In this calculation, it is useful to realize that u^d^ = d T , i.e. the proper 
time derivative. Then, we calculate that c* M <9 r £ = d T d^. Furthermore, for 
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adiabatic changes, cf. Eqs. (^J), we make use of the relation, 



d (e + p\ = dp 
\ n J n 

With these ingredients, we finally obtain: 

(e + + u ft u u d l/ p = d^p , (14) 

which is the hydrodynamic four-vector equation of motion. 

Separating the space and time components of Eq. (|l4|), recalling Eqs. (|J), 
we immediately obtain the results: 

[dt+v-Vh 

[dt + v • V]t7 

i.e. the relativistic energy and momentum flow equations, the latter being 
known as the SR generalization of the Euler equationJ13 



1 



1 



(dtp-j^dt+v- V]p) 



7 2 (e +P) 



(Vp + vd t p) 



(15) 
(16) 



2.3 Example: Relativistic Rayleigh-Plesset Equation for a Spherical Bubble 

In all the applications of relativistic hydrodynamics to multiparticle produce 
tion relatively little attention has been payed to the spherical geometry.tMlia 
The longitudinal-expansion is usually considered predominant compared to the 
transverse one.Bt-3. However, in collisions of very heavy nuclei the transverse 
expansion might be equally important. Therefore, applying the variational 
principle, we derive an effective equation of motion for a spherically symmetric 
fluid bubble. We sketch the most simplified scenario, while more details can 
be found elsewhere Jld 

Let no(r,t) denote the baryon number density inside a bubble of radius 
R(t) for a space fixed coordinate system (center-of-mass) . In a comoving frame 
the density is n = 7~ 1 no, with 7 _1 = (1 — v 2 ) 1 / 2 . We make the ansatz, 

n (r,t) = R- 3 f(x- ia (t)), x = r/R , (17) 

where / denotes the shape function depending parametrically on a (or several 
components a^) and on the scaling variable x. Then, the continuity equation, 
dtriQ + r~ 2 d r (r 2 vno) — 0, can be solved for the radial velocity, 

v(x,t) = xR - ^ [ dx' x' 2 ^J- . (18) 
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Since the constraints incorporated in the action Im , Eq. (g) , are automatically 
satisfied by now, the effective action becomes: 



l RP 



dt L = —Am dt R 3 dx x 2 e (f(x; a)/( 7 i? 3 )) . (19) 
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The corresponding Euler-Lagrange equations present coupled ordinary differ- 
ential equations for R(t) and a(t), which demonstrates the essential simplifi- 
cation due to the variational ansatz. 

In the limiting case of homologous motion, i.e. a time independent shape 
function f{x) and consequently a linear velocity profile, v(x,t) = xR, the 
resulting Rayleigh-Plesset type equation is: 

^j^R 3 J dxx A 1 2 {p + e)^=m 2 dxx 2 p , (20) 

which fully incorporates the necessary SR corrections. 

Assuming an ultrarelativistic adiabatic equation of state for an idealized 
quark-gluon plasma fluid, we obtain the relations, 

p = Po [R /R(t)} 3r f r (x) 1 - r , £ = s p/p , (21) 

with the adiabatic index T = 4/3. The nonperturbative QCDjj/acuum pressure 
B w (200MeV) 4 acting on such a plasma from the outside,0 is incorporated 
here by replacing p — > p — B on the right hand side of Eq. (20). Then, it is 
straightforward to solve Eq. (|2(]) together with Eq_(E[l]) numerically for suitable 
initial conditions, which illustrates our approach.!!!! 

Along these lines, also an expanding hadronic phase outside the plasma 
bubble and dissipative processes due to the phase boundary can be incorpo- 
rated and are presently under study. 



3 General Relativistic Hydrodynamics 

For most astrophysical or cosmological applications of the variational approach 
it will be necessary to incorporate GR gravity. According to our considerations 
in Sec. 2.2, we have to combine the usual action for gravity, Iq, with a suitable 
fluid/matter action, Im- The following total action suggests itself: 

I G + I M = J d^x^f—g | ~e(n)+ £(x) V M [n«"] + ~£(x) - 1] J , 

(22) 

where the only difference in the matter part, cf. Eq. (B), arises from co- 
/invariance requirements; we use standard notation.l!3 Here the variational 
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principle is: SS — 0, for all independent variations of the metric g^ and of 
m m , n, £, C, as before. 

Performing the variations, using 9 M In y 7 — 5 = r^ Q : we obtain Eqs. @, 
(|H), and (H). Furthermore, varying £ and g^ v leads to, respectively: 

v,K]E(9^r;)K] = o, (23) 

(SnGy'G^ - i^Zp^+p)^ = , (24) 

where is the Einstein tensor, and we used d g u In \J—g = —g^/2. 

Identifying the energy-momentum tensor by Slu/^g^v = \J — /2, we 
observe that Eq. ( p4| ) is nothing but Einstein's equation for an ideal fluid, 
G* v = 87rGT^, with = (e + p)u^u v - pg^ . Then, the Bianchi identity, 
V M G^ = 0, implies the GR hydrodynamic equations: V = OH 

We conclude by briefly presenting the related effective Lagrangian for the 
most general spherically symmetric system, where the metric is: 

ds 2 = c^dt 2 - e 2X dr 2 - R 2 dn 2 , (25) 

with cf>, A, R unknown (variational) functions of r, t. Choosing the comoving 
frame with it M = (u°, 0, 0, 0), the normalization u 2 = 1 implying u° = exp(— <fr), 
and with the (baryon) density n being proportional to i? _2 exp(— A), the cur- 
rent conservation follows. Thus, the constraints being satisfied, the Lagrangian 



corresponding to Eq. (22) for spherical symmetry is obtained: 



£ = e" i c A i? 2 [-e(n) + \—. ,{R'[2R4>' +R']e- 2X -R\2R\+R]e~ 2,i > + \}] , (26) 

\\S-n(jR l 

where / = df/dt and /' = df/dr; we omitted a total derivative which does 
not influence the related Euler-Lagrange equations for (j>, A, R. 

We have demonstrated that known cases like the Misner-Sharp equa- 
tion relevant for stellar collapse or explosion, the Tolman-Oppenheimer-Volkov 
equation describing neutron star structure, and the Friedmann-Robertson- 
Walker cosmological model can be re-derived from this action. 



4 Conclusions 

All known results involving ideal fluids in SR and GR, respectively, can be 
obtained directly from the actions, Eqs. (ph and (p2|) , respectively which we 
derived. Some earlier attempts can be found in the literature.E£l However, 
we motivated in Sec. 1 and stress here once more that the application of the 
variational principle together with physically motivated trial functions, as il- 
lustrated in Sec. 2.3, may provide a powerful tool to study new and interesting 
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fluids under extreme conditions, such as in ultrarelativistic nuclear collisions, 
or some more complex systems in GR. We will report the details of our calcu- 
lations and such new applications elsewhere. 
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